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Abstract 

The generalized spherical Radon transform associates the mean val- 
ues over spherical tori to a function / defined on § 3 C H, where the 
elements of S 3 are considered as quaternions representing rotations. It 
is introduced into the analysis of crystallographic preferred orientation 
and identified with the probability density function corresponding to the 
angle distribution function W. Eventually, this communication suggests 
a new approach to recover an approximation of / from data sampling W. 
At the same time it provides additional clarification of a recently sug- 
gested method applying reproducing kernels and radial basis functions 
by instructive insight in its involved geometry. The focus is on the cor- 
respondence of geometrical and group features but not on the mapping 
of functions and their spaces. 

1 Motivation and Introduction 

In texture analysis, i.e. the analysis of preferred crystallographic orien- 
tation, the orientation probability density function / representing the 
probability law of random orientations of crystal grains by volume is a 
major issue. In x-ray or neutron diffraction experiments spherical in- 
tensity distributions are measured which can be interpreted in terms 
of spherical probability distributions of distinghuished crystallographic 
axes. In texture analysis, they are referred to as pole probability density 
functions. Generally, if / is the orientation probability density function 
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of a random rotation represented by its random quaternion variable Q, 
then the probability density function of the random direction QhQ*, 
where h £ S 2 C R 3 is a fixed crystallographic direction, is provided by 
the 1 -dimensional spherical Radon transform, integrating / over all 1- 
dimensional great circles C C S 3 representing all rotations mapping h 
onto a direction r £ S 2 . 

Extracting pole densities from raw intensity data collected in diffrac- 
tion experiments requires several steps of processing, in particular defin- 
ing the "background intensity" and correspondingly normalizing the in- 
tensity data. Any step in this process is largely subject to primary ex- 
perimental errors. The influence of experimental errors can easily be 
reduced by considering mean intensities resulting from an additional 
rotation of the specimen during the measurements (Nikolayev, 2004). 
These mean values correspond to integration of spherical pole density 
functions over small circles. Thus, they correspond to mean values of 
/ over 2-dimensional spherical tori T C § 3 with core C (Meister and 
Schaeben, 2004) which are called generalized spherical Radon transform 
(Helgason, 1994; 1999). Here, it is identified with the angle density func- 
tion W (cf. Bunge, 1969, 1982). Independently of the new experimental 
possibilities, they have proven to be instrumental for the inverse spheri- 
cal Radon transform (Muller et al., 1981; Helgason, 1994; 1999). 

2 Geometry of Rotations 

One of the most beautiful features of quaternions is the role they play in 
the representation of the rotations of the low dimensional spaces M 3 and 
M 4 . The description here is based on the book by Delanghe, Sommen 
and Soucek (1992). 

2.1 Rotations and Quaternions 

The skew-field of quaternions EI is generated by the elements 
ej, i — 1, 2, 3, which fulfil the relations 

(i) e 2 = -1, i = 1, 2, 3; 

(ii) eie 2 = e 3 , e 2 e 3 = e ± , e 3 ei = e 2 ; 

(iii) titj + ejSi = 0, i,j = 1, 2, 3; i ^ j. 
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The unit element of EI should be denoted by 6q. Then an arbitrary 
quaternion q £ EI can be represented as 

3 

q = g + q = <?oe + ^ q^i = Sc q + Vec q, 

i=i 

where q = Scq is the scalar part and q = Vecg is the vector part of 
q. If Scg = 0, then q is called a pure quaternion, the subset of all pure 
quaternions is denoted by Vec EI. The subset of all quaternions with 
vanishing vector part may be denoted by Sc EI. In this way R 3 and R are 
embedded in H and R 4 = R 3 © R may be identified with EI. 

The quaternion 

3 

Q* = Qo - q = <?oe - ^2 Qiei = Sc q - Vec q 

i=l 

is called the conjugate of q. The conjugation is an anti-involution, i.e. 

(pq)* = q*p*, and 

3 

* * II 1 1 2 \ ^ 2 

qq = q q = \\q\\ = ^Qj, 

j=0 

where \ \q\ \ is the norm of the quaternion q which coincides with the norm 
of q in the Euclidean space M 4 . We denote the set of all unit quaternions, 
i.e. all quaternions with norm 1, by § 3 . 

A linear transform T in M 3 or R 4 is called a rotation if and only if it 
leaves the scalar product in the Euclidean space and the orientation 
invariant. Of course we have that T is invertible with detT = 1. The 
set of all rotations (in the previous sense) in R 3 and R 4 is given by the 
special orthogonal group SO (3) and SO (4) respectively. 

Proposition 1. Let q £ S 3 and define the linear transform T q : R 3 — > 
R 3 by 

T q (x) = qxq*. 

Then T £ SO(3) if and only if T = T q for some q £ § 3 . 

Thus any rotation in R 3 can be represented by a unit quaternion q via 
qxq* . As it turns out, quaternions can also be used to represent rotations 
in R 4 , i.e. the group SO(4). 
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Proposition 2. Let q, p G S 3 be unit quaternions and define the maps 

$ g , \&p : M 4 — > M 4 by = gx and ^ p {x) = xp. 

Then T G SO (4) if and only if there exist q, p G § 3 such that T = Qq )P *, 
where 

Thus any rotation T G SO(4) may be represented by T(x) = qxp*, 
where q, p G S 3 . 

Proposition 3. Let be SO (4) 9 T = 9 3)9 » = T 9 . Then T G SO (3). 
Obviously, 

T(e ) = O g ,g*(e ) = r g (eo) = qe q* = qq* = 1 = e , 

i.e. T fixes the North Pole eo and hence T G SO (3). 

Remark 1. It is easily seen that § 3 is a double covering o/SO(3) and 
S 3 x S 3 is a double covering o/SO(4). 

Representing rotations by quaternions yields an instructive and geo- 
metrically appealing clarification of the geometry of rotations. Here we 
summarize the major results of Meister and Schaeben (2004) and add 
some explicit relationships which will prove helpful for our purposes. 

2.2 Geometrical objects of § 3 C i and of § 2 C M 3 

Definition 1. Let q\,qi G S 3 be two orthogonal unit quaternions. The 
set of quaternions 

q(t) = q 1 cost + q 2 smt, £e[0,27r), 
constitutes the great-circle denoted C(qi,q2) C S 3 . 
With 

1 — rh h + r , s 

qi = 77- 7-77 and q 2 



|1 — rh\\ \\h + r\\ 

for h, r G § 2 with r ^ — h it is obviously 

bc(q{t)) = cos — — = cos — cos t. 

Obviously, (h, r) and (— h, — r) define the same great circle 
C(9i,? 2 )c§ 3 . 
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Definition 2. Let qx,q2,Q3,Q4 £ §> 3 be four mutually orthonormal 
quaternions; let C(q\, (72) denote the circle spanned by quaternions qi,q2, 
and C(q3,q4) the orthogonal circle spanned by q3,q4- The set of quater- 
nions 

q(s,t;Q) = (qi cos s + q2 sin s j cos 6 + (q% cos t + q^ sin t j sin , 

s,t G [0,2tt), 6 G [0,tt/2] (2) 

constitutes the spherical torus denoted T(C(q\, 92); 0) C S 3 u^/j core 

C(?l,?2)- 

Definition 3. Lei r G § 2 be a unit vector and 

r(t) = cos — h r sin -, t G [0, 27r), 
represent the rotation about vbytE [0, 2ir). Then the set of unit vectors 

r f (t) = r(t)T f r*(t), tG[0,2vr), (3) 

with r' G § 2 in the plane spanned by h and r snc/i i/iai r -r' = cos(p), h • 
Tq = cos(?7 — p) constitutes the small circle or cone 

c(r; p) = {r' G § 2 | r • r' = cosp} C § 2 

wt/i ang/e p itrat/i respect to its centre r. 

Its parametrized form explicitly reads (e.g. Altmann, 1986) 

r'(t) = r' cost + (r x r' ) sint+ (r • r' ) r (1 - cost), t G [0,2tt). 

With h G S 2 and 

h(t) = cos^ + hsin^, tG[0,27r), 

representing the rotation about h by t G [0, 27r) and with hp G § 2 in the 
plane spanned by h and r such that h • = cos(p), r ■ h' Q = cos(r/ + p) 
analogously results in 

h'(t) = h(t)h(,/i*(t), t G [0,2tt). (4) 
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2.3 Fibres 

The fibre G(h, r) C SO (3) of all rotations with = r is represented by 
the circle C(qi,q 2 ) C § 3 spanned by unit quaternions q±,q 2 G S 3 given 
in terms of h, r G § 2 by Eqs. (JTJ), for example. Therefore, the notation 
Ch, r = C(qi,q 2 ) is used where it is more instructive keeping in mind 
that Ch r = C_h,- r - Thus the major property of the circle C(qi,q 2 ) is 
that it consists of all quaternions q(t),t G [0, 2ir), with q(t)hq*(t) = r 
for all t G [0, 2w), and that it covers the fibre G(h, r) twice; moreover, it 
is uniquely characterized by the pair (h, r) G § 2 x §> 2 and its antipodally 
symmetric (— h, — r). 

The set of all rotations mapping h on the small circle c(r; p) is equal 
to the set of all rotations mapping all elements of the small circle c(h; p) 
onto r and represented by the spherical torus T(C(qi, q 2 ); f) C S 3 with 
core C(qi,q 2 ). 

The distance d of an arbitrary q G § 3 from the circle C(q\, q 2 ) is given 

by 

d(q, C(qi, qa)) = - arccos(ghg* • r) 

If d(q, C(qi, q 2 ) = p, then q and C are called p-incident. 

Then, the torus T(C(qi, q 2 ); |) consisting of all quaternions with dis- 
tance | from C(qi,q 2 ) essentially consists of all circles with distance | 
from C(qi, q 2 ) representing all rotations mapping h on c(r; p) and map- 
ping c(h; p) on r. 

The representation of a torus can be factorized in the following way. 
Let (h, u), (r, v) G § 2 x § 2 with u ■ v = cos p. Then G x = C h ,u C § 3 and 
C 2 = C r v C § 3 exist such that 

Pihp* = ri for all pi G C\ 
p 2 h 2 p* 2 = r 2 for all p 2 G C 2 

Then p 2 p\ maps h on the small circle with centre r and angle p = 
arccos(u • v) as 

P2Pihp{p 2 ■ r = pihpi ■ p 2 rp* 2 = u • v, (5) 

i.e. {p 2 pi \pi G Cj, i — 1,2} represents the torus T(C;p/2) with core 
C = C*h,r and angle p = arccos(u-v). Letp G T(Ch, r ; p/2), thenphp* = x 
with x ■ r = cos p. Given u, v G S 2 with u ■ v = cos p there exists a unique 
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P2 G S 3 such that P2W2 = u an d Vi^p\ = v - For any p\ G Ch, u it is 
p\p2php*p2P\ = h implying p = p\p\ . 

Since g(t)hg*(t) = r implies that q(t)h'q*(t) • r = h' ■ q*(t)rq(t) = 
h'h for all q{t) G C(qi,q 2 ), the set {q(t)h'q*(t)} represents the small 
circle (cone) around r with angle t = arccoshh'. In parametrized form 
employing Eqs. (jHJ), (0J) 

q(t) h'(u) q*(t) = r'(u + 2t), t,u G [0, 2vr), q(t) G C( 9l , tfe) (6) 

(Meister and Schaeben, 2005), which may be rewritten for quaternions 
as q(t) = r(u + 2t). 

Then the distance of q(t) G C(q±, q^) from an arbitrary circle C\ repre- 
senting all rotations mapping h x on V\ is given by spherical trigonometry 

d{q{t),d) = - axccos(g(t)hig*(t) • n) (7) 

= - arccos ^hhi + rri + y/ (1 — (hhi) 2 )(l — (rri) 2 ) cos . 

Eventually, the set of all circles C(pi,p2) C S 3 with a fixed distance 
% of a given q G S 3 , i.e. the set of all circles tangential to the sphere 
s(q; p/2) with centre q and radius p/2, is characterized by 

I = d(<7,C(pi,p 2 )) = - arccos (ghg* • rj, (8) 

where r G § 2 is uniquely defined in terms of h and pi,P2 by r := 
p(t)hp*(t) for all p(t) G C(p±,p2) and any arbitrary h G § 2 , i.e. each 
circle C(p±,p2) represents all rotations mapping some h G § 2 onto an 
element of the small circle c(qhq*; p). Thus, for each q G § 3 and p G [0, 7r) 

{c7( Pl ,ft)|d( g ,C7(p 1) p 2 )) = |}= |J |J C h , r (9) 

heS 2 r£c(qhq*;p) 

= U U <?(Pi(h,r) I p a (h,r)) 

hGS 2 rgc(<j h <jr* ;p) 

= U U C(Pi(h,r),p 2 (h,r)), 

hGS^ r€c(<jhq*;p) 

where denotes the upper hemisphere of S> 2 . The last equation is due 
to the fact that (h, r) and (— h, — r) characterize the same great circle 

C*h,r = C-h,-r- 
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3 Radon transforms 



3.1 The spherical Radon and generalized spherical Radon 
transform 

Experimentally accessible crystallographic pole figures P(h, r) of 
recorded diffracted x-ray or neutron beam intensities are modelled by 

P(h.r) = -L f f(g)dg, (10) 

47r </G(h,r)UG(-h,r) 

and interpreted as the probability density function that the crystallo- 
graphic axis ±h statistically coincides with the specimen direction r 
given the orientation probability density function / : SO (3) h- > 

Let C denote the set of all 1-dimensional totally geodesic submani- 
folds C C § 3 . Each C £ C is a 1-sphere, i.e. a circle with centre O. Each 
circle is characterized by a unique pair of unit vectors (h, r) £ § 2 x § 2 
(and equivalently by its antipodally symmetric) by virtue of q h q* = r 
for all q £ C. Thus, refering to the quaternion representation, Eq. (|10j) 
can be rewritten in a parametric form as 

P(h,r) = -±- f f(q)duj 1 (q), 

where the circle CcS 3 represents all rotations mapping h £ S 2 on r £ 
S 2 and where C (qi, q-z) is the orthogonal circle representing all rotations 
mapping — h on r, and where U\ denotes the usual one-dimensional 
circular Riemann measure. 

Following Helgason (1994; 1999), 

Definition 4. 

±- [ f(q)duj 1 (q)= [ f(q)dm(q) = (Kf)(C) 
^ Jc Jc 

with the normalized measure m = is referred to as 1-dimensional 
spherical (totally geodesic) Radon transform of f whenever f is integrable 
on each great circle. 

The Radon transform of / may be represented as the convolution of / 
with the indicator function of the great circle C. It associates with the 
function / its mean values over great circles C £ C. Since each great 



9 



circle is uniquely characterized by a pair (h, r) £ S 2 x § 2 and its antipo- 
dally symmetric, we also use the notation (7£/)(h, r) = (1Zf)(— h, — r) 
whenever it is more instructive. 

It should be noted that no distinction has been made whether / refers 
to SO(3) or § 3 , even though the form of / depends on the representation 
of g\ in particular, with respect to § 3 only even functions / could be 
orientation probability density functions as q £ § 3 and — q represent the 
same orientation. Then 

P(h,r) = ^((ft/)(C h , r ) + (ft/)(C_ h ,r)) = W)(h,r), (11) 

where X f is referred to as basic crystallographic x-ray transform. 

Further following Helgason (1994; 1999) the generalized 1- 
dimensional spherical Radon transform and the respective dual is well 
defined. 

Definition 5. The generalized 1 -dimensional spherical Radon trans- 
form of a real function f : S 3 > M 1 is defined as 

(n^f)(C) = -^— [ f(q)dq. 
4vHsmp J d ( q ,c)=p 

The generalized Radon transform of / may be represented as the con- 
volution of / with the indicator function of the torus T. It associates 
with / its mean values over the torus T(C, p) with core C and radius p, 
Eq. ©. 

Following Berens et al. (1968), Freeden et al. (1998, p. 64), 

Definition 6. The spherically generalized translation of a real function 
F £ C(§ 2 ) or F £ L p (S 2 ), 1 < p < oo, is defined as 

(TV>F)(r) = 1 / F(r')dr'. 

27Ta/1 — COS^p Ar'=cosp 

where 1tx^/\ — cos 2 p is the length of the circle c(r; p) centered at r with 
radius cos p. 

It can be determined by 

(T^F)(r) = — / F(v'(t))dt, 

I'K-J 1 — COS^ p Jo 
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with r'(t) given according to Eq. Q. 

When the translation is applied to the Radon transform with re- 
spect to one of its arguments, then the geometry of rotations represented 
by quaternions amounts to 

(TM[fc/]Wr) = -i— / (Kf)(ti,v)dh' 
\ J Iirsmp J c{h . p) 

= 77^— [ W)(h,r')dr' (12) 
2nsmp J c(r . p) 



1 



c(r;p) 

f{q)du 1 {q)dr' 



Alt 2 Sinp i c ( r;p ) Jc( qi (h,r'),q 2 (h,r')) 

^— I f(q)dq (13) 

' SUl p J T ( C (a, (h,rWh,r));-£) 



AlT 2 Smp JT(C(qi(h,r),q 2 (h,r)y^) 

All 2 smp J d( g> c(gi(h,r), g2 (h,r)))= 
(^^ /2) /)(C , h ,r). (14) 



f(q)dq 



Eq. ((EI), cf. (Bunge, 1969, p. 47; Bunge, 1982, p. 76), is an Asgeirsson- 
type mean value theorem (cf. Asgeirsson, 1937; John, 1938) justifying 
the application of to IZf regardless of the order of its arguments, 
and Eq. (jlHjl is instrumental to the inversion of the spherical Radon 
transform Helgason, 1994; 1999). We have just accomplished 

Proposition 4. The generalized 1- dimensional spherical Radon trans- 
form is equal to the translated spherical Radon transform 

V^[^/])(h,r) = (^ 2 V)(C h ,r) 
and it can be identified with the angle density function 

W)(h,r;p):=-^— / (Kf)(h, v')dx' (15) 
27rsmp J c{r . p) 

The angle density function (Af)(h, r; p) has been introduced in (Bunge, 
1969, p. 44; Bunge, 1982, p. 74) (with a false normalization). It is the 
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probability density that the crystallographic direction h statistically en- 
closes the angle p, < p < n, with the specimen direction r given the 
orientation probability density function /. It should be noted that 

(Af)(h,r;0) = (ft/)(h,r), (16) 
W)(h,r;7r) = (Kf)(h,-r). 

Finally, with respect to the diffraction experiment, it should be no- 
ticed that 

(T P [X /]) (h, p) = i ((n^f)(C) + (TZ^f)^)) 

\((Af)(h, r; p) + (Af)(-h, r; p)) 
= (W)(h, r;p) (17) 

Thus with respect to a diffraction experiment, (W)(h, r; p) is at once 
accessible if the specimen rapidly rotates around the specimen direction 
r during the measurements. 



3.2 Kernels and their twofold Radon transform 

Now let K be a kernel function S 3 x S 3 i— > Ri. Then we may apply the 
Radon transform twice in the sense that we apply it once with respect 
to the first and once with respect to the second variable (cf. Boogaart 
et al, 2005), i.e. 

(KiK^p,)})^) = i- I K(p 1 ,p 2 )du 1 (p 1 ) = F(C 1 ,p 2 ) 

and 

(nKKiK^o,)})^)})^) = {n[F{c u o 2 )]){c,) 

= t~9 / K(p l ,p 2 )du) 1 {p l )du l (p 2 ) = G(Ci,C 2 ), 
4tt 2 Jc 2 Jc! 

where C\ denotes the great circle representing all rotations mapping h x 
on Pi, and analogously C 2 with respect to h 2 and r 2 . 
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In case the kernel K depends only on the quaternion product p\p2, 
then 

(^[(^[iTKo 2 )])(C 1 )])(C 2 ) = 7^/ / Kiplp^dooMdooM 

K{p) dp 



4vT 2 a/1 - (ri ■ r 2 ) 2 ./(phip*-h 2 )=(ri-r 2 ) 

4vr 2 sm(p/2) i T ( C ;p/2) 
= (^ ( " /2) JO(C) 
= (AK)(h u h 2 ;p) 

with the core C = Ch 1 ,h 2 an d P — arccos(ri • r 2 ) according to Eq. (J5J). 

In case the kernel K depends only on uj/2 = arccos(Sc(p^p 2 )), then 
the function F depends only on arccos(p 2 hip 2 -ri) = d(p 2 , C\). As a func- 
tion of the variable p 2 with parameters hi, implicitly provided by Ci, 
F may be referred to as fibre function. They share most of the charac- 
teristics of ridge functions initially defined in linear spaces and discussed 
in (Donoho, 2000). Then the function G depends only on arccos(h!h 2 ) 
and arccos(rir 2 ), cf. (Boogaart et al., 2005). 

Proposition 5. For a kernel function of the form K(pi,p 2 ) = K(p\p2), 
its angle density function, i.e. its generalized spherical Radon transform, 
is identical with the twofold application of the spherical Radon transform 
with respect to the two components of its argument. 



3.3 The dual spherical Radon and dual generalized spherical 
Radon transform 

Following again Helgason (1994, 1999) we have two more definitions. 

Definition 7. The dual l-dimensional spherical Radon transform of a 
real continuous function </? : C i— > R 1 is defined as 

JC3q 

where \i denotes the unique measure on the compact space C G C : q G C , 
invariant under all rotations around q, and having total measure 1 . 
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Thus (IZcp) (q) is the mean value of ip over the set of circles C G C passing 
through q. 

Definition 8. The dual generalized l-dimensional spherical Radon 
transform of a real continuous function ip : C i— ► M 1 is defined as 



(fc<"V)(?) 



^(C)dJl(C). 



{C£C:d{q,C)=p} 



With the normalized measure Ji = A~ 1 (p)dfi, A(p) = An sin 2 p, it is the 
mean value of tp over the set of all circles C G C with distance p from q. 
Again, geometrical reasoning, Eqs. (JBJ), ©, yields that 

\TZ^ (1Zf)(o) j(q) is the mean value of (11 f) over the set of all 

circles C G C with distance p from q, i.e. tangential to the sphere s(q; p). 
More specifically, with the usual two-dimensional spherical Riemann 
measure u>2, 



n { P ) 



(llf)(o) )(q)= / (lZf)(C h , r )d^(C h 

d(q,C)=p 
1 



47rsin(p/2) J S 2 J c{qhq * ]p/ 2) 



(llf)(C Kr )drduj 2 (h) 



1 

2 ./ 

1 

2 



(Af)(h,qhq*;p/2)du; 2 (h), 
(K^f)(C h , qhq * du 2 (h) 



K 



W)(o) (q) 



(18) 
(19) 

(20) 
(21) 



which is instrumental for the inversion of the spherical Radon transforms. 
In particular, 

(£[(ft/)(o)])(g) = 1 -jjlZf)(h,qhq*)duj 2 (h), (22) 

Next we want to describe the generalized dual Radon transform in 
a more group-theoretical way (cf. Helgason, 1994; 1999). Let Co G C be 
a fixed great circle with distance p from q, i.e. d(q, Co) = p. Let Q r)S * 
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be the rotation that maps the North Pole eo to q, i.e. q = reos* = rs*. 
Since the distance on S 3 is rotational invariant we obtain 



Since all rotations in SO (3) leave the North Pole eo invariant we conclude 

d(q, Co) = d(eo, r*Cos) = d(keok*, kr*Cosk*) = d(eo, kr*Cosk*), 

hence the set {kr*C sk*, r, s G S 3 fixed, k G § 3 } is the set of all great 
circles that have distance p from the North Pole e$. Now, the North Pole 
eo is rotated back into q leading to 

d(q, Co) = d(e , kr*C sk*) = d(re s*, rkr*C sk*s*) = d(q, rkr*C sk*s*). 
Thus the set of all great circles with distance p from q is given by 



The previous consideration amounts to a new description of the gen- 
eralized dual Radon transform: 



where dk = ^ sin (3 da d/3 d'j is the invariant Haar measure on SO (3), 
Co G C a fixed great circle such that d(q, Co) = p and q = rs*. 

Analogously to the spherically generalized translation for function on 
S 2 , we have 

Definition 9. The spherically generalized translation of a real function 
f G C(S 3 ) or F G L P (S 3 ),1 < p < oo, is defined by the mean value 
operator 



where A(p) denotes the surface area of the sphere s(q; p) centered at 
q G S 3 with radius cos p. 

Then the following proposition can be shown. 



d(q, C ) = d(rs*, C ) = d(e , r*C s). 



{rkr*C sk*s*, k G § 3 ,g = rs*}. 




SO (3) 
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Proposition 6. 



W)(o) )(q) 



(Kf)(C) dfl(C) 



{C&C:d(q,C)=p} 

(T d ^f){q) dm(p), 



(23) 



Co 



where Cq is a fixed great circle with distance d(q, Co) = p from q. 



W)(o) )(q) 



(Kf)(C) dfl(C) 

{C£C:d(q,C)=p} 

(Kf)(rkr*C sk*s*)dk 



SO(3) 



SO(3) JCo 



f(rkr*psk*s*) dm(p) dk 



f(rkr*psk*s*) dk dm(p) 

C JSO(3) 

While k varies in SO(3), the set {rkrpsk* s* , k e SO(3)} is a small 
sphere around q containing p : 

d(q, rkrpsk* s*) = d(rs*, rkrpsk* s*) = d(eo, krpsk*) 
= d(k*e k, r*ps) = d(e , r*ps) = d(rs*, p) = d(q, p) , 

and dk = A~ 1 (p)dp is the normalized Riemann measure of the small 
sphere s(q; p) with center q and radius p = d(q,p). Hence 



n { P ) 



W)(o) )(q) 



Co JSO(3) 
1 



Co 



f(rkr*psk*s*) dk dm(p) 
f(p) dpdm(p) 



S P (q) 



(T d ^f)(q)dm(p). 



(24) 



Co 



Remark 2. A special situation occurs when f is a central function. 
Assuming that f is central with respect to q , (TZf) is constant on the 
set of all great circles C G C for which q Q h. c qQ • Tc> where h c , r c denote 
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the unit vectors associated to the great circle C. This is true for the set 
of all great circles with d(qo, C) = p as they are characterized by Eq. (0). 
Thus we can drop the outer integration and find for the right hand side 
of Eg. 



W)(o) (g ) = W){C) for some C with d(q , C) = p. 



n ( P ) 



Moreover, in this case, {T p f){qo) = /((Zo) for all p. Thus the left hand 
side of Eq. 



(T d ^f)(q ) dm(p) = (Kf)(C) for some C with d{q ,C) = p. 

d(qo,C)=p 

Parameterizing the great circle Cq by polar coordinates and assume 
/ even (f(q) = /(—<?)) Eq- (EH) takes the form 



ftp) 



W)(o) )(q) 



(T d ( q 'rif){q) dr 



(25) 



Denote by q G C a point with minimum distance from q, i.e. d(q , q) = 
p. Using spherical trigonometry with respect to the triangle qqop we 
obtain 

cos d(q,p) = cosd(qo,q) cosd(qo, p). 
Fix q and set v = cosd(qo, q) and u = v cosd(qo, p) = cos d(q,p), 



n (p) 



F(u) = (T d ^f)(q), F(v) 

Then Eq. (J2*5*j) becomes an Abel's integral equation 

F(u 



W)(o) )(q). 



F(v) 



which is inverted by 
F(u) 



1 d 

2tt du 



o \/v 2 — u 2 



F(v) 



du, 



dv. 



Since F(l) = lim u _^ 1 (T arccos "/)(g) = f(q), we get the following inversion 
formula (cf. Helgason (1999)) 



'<»> = T, 



d I I ^(arccosti) 

du 



W)(o) )(q) 



\Ju 2 — V 2 



dv 



u=l 
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and with Eq. (|19|) we obtain 



■7- / / (Af) (h, qhq*; 2 arccos v) du 2 (h) 7= = dv 
du J J S 2 Vn 2 - v 2 



u=l 

which will be further transformed. First, we substitute t = v 2 , getting 

1 d r 2 f , . n 1 



^7T du 



[ [ (Af) (h, qhq* ; 2 arccos y/t) dm 2 (h) dt 

Jo Js 2 y/u 2 — t 

1 d f u2 f 1 
~a i i (Af)(h, qhq*; 2 arccos y/t) du 2 (h) - dt 

47T du 2 J J S 2 y/u 2 - t 



u=l 



u=l 



Next, we put s = u to get 



= - — f [ (Af) (h, qhq* ; 2 arccos y/i) du 2 (h) — ; dt 

47T ds J Q J S 2 y/s-t 

To shift the singularity, we set w = s — t which gives 

/ / (Af)(h, qhq*; 2 arccos y/s — w) duj 2 (h)—= dw 
Jo i§2 y/w 



in ds 

Now, differentiation gives 



1 

47T 



(Af)(h,qhq*;7r) du 2 (h) 

[ f -^-(Af)Ch, qhq*; 2 arccos y/s — w) dh—= 
Jo is 2 ds V 1 



s = l 



Using -£(Af) = --^{Af) and taking into account (Af)(h, qhq*; n) 
(Kf)(h, -qhq*) (cf. Eq. and s = 1 result in 

= ^(JW)(h,-qhq*)dou 2 (h) 

~ f f ~r~ ((Af)(h,qhq*; 2 arccos y/T^w)) du 2 (h)-^=dw 

J J S 2 dlU y/W 

and finally put 2w = 1 — cos 9 to obtain 
(Kf)(h,-qhq*)dcu 2 (h) 

+ 2 r I -j^-A{Af){h )q hq*;9)) dw 2 {h)cos\dd] 
J J S 2 dcos6 ) 



1 

Air 
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4 Inversion of the generalized spherical Radon 
transform 

For / £ C°°(§ 3 ) the spherical Radon transformation 1Z : f \— > IZf has 
a kernel consisting of the odd functions, i.e. the functions satisfying 
f(q) + f(—q) = 0. Since q and — q represent the same rotation, and since 
/ is a probability density function, we are interested in even functions 
/ : S 3 i— > IR+ only. In practical texture analysis, we are especially inter- 
ested to recover / from data u>j = iw(hj, r^;pj) originating in sampling 
(W/)(h, r;p) for discrete values of h, r £ S 2 and p £ (0, 7r/2). Obvi- 
ously, due to the additional symmetry Eqs. (fTTj). ??), introduced by 
Friedel's law, / cannot completely be recovered from integral intensity 
data Li = P{hi, r$) nor by mean integral intensity data Wi = w(hi, r«; pi). 
In terms of spherical harmonics, only harmonic coeffficients with respect 
to harmonics of even order can be determined from diffraction intensity 
data (Matthies, 1979). 

4.1 Inversion of the spherical Radon transform 

In texture analysis, i.e. in material science and engineering, the best 
known inversion formula dates back right to the begining of "quantita- 
tive" texture anaylsis (Bunge, 1965; Roe, 1965). The formula may be 
rewritten in a rather abstract way as 

/ = ^so{3) $ -T^xs 2 Tlf , (26) 

where S denotes a scaling matrix with entries \/2£ + 1 indicating the ill— 
posedness of the inverse problem. It states that the harmonic coefficients 
of / are up to a scaling equal to the harmonic coefficiens of its Radon 
transform. In the context of texture analysis and experimentally acces- 
sible "pole figures", the former statement is true only for even-order 
coefficients. 

The first analytical inversion formula was conributed by Matthies 
(1979) and then rewritten in terms of the angle density fucntion (Muller 
et al, 1981). 
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In integral geometry the following inversion formlae are known 

/ = J- / (-2A s2xS2 + 1) 1/2 (ft/)(h, gh) duj 2 (h) 
Air J S 2 

= ^-jz[(-2A s , xS2 + lf/ 2 (TZf)} (27) 

where 71 denotes the adjoint operator of TZ with respect to L 2 , and with 
respect to Eq. (J22J) 

/ = J_(-4Ass + i) 1/2 £rc/ 

= i-^Aga + l) 1 / 2 ^)^/ (28) 
= ^(-4A s a + If/ 2 (nKf) (29) 

cf. (Helgason, 1994; 1999). The equivalence of all these formula is shown 
in (Bernstein and Schaeben, 2005). 

Comparing Eq. (|2*5j) with Eq. (}2"9~j) with respect to Eq. (j2T?j) we may 

ask ourselves: Is ^ 

(-4A S 3 + I) 1 ' 2 (KKf) 

a better conditioned inverse problem in practical applications than 

((^Aga + l) 1 / 2 ^) TZf 
which is conditioned like 

^ n F SO(3) ^ F 52 x§2 

Obviously, (— 4AS3 + 1) 1 / 2 is worse conditioned than (— 4A§3 + l) 1//2 7?., and 
1Z1Z is better conditioned than 7Z. Thus, the question can be specified 
whether the differences just cancel out, or does the later improvement 
dominate? This problem will be pursued elsewhere. 
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4.2 Harmonic series expansion 

In terms of spherical hamonics, Eqs. (??, I15JI result in 

W)(h,r;p) = 1 / ^^Q m T(h)F;(rVr' 

2^1- cos 2 p J c{r;p) , =2/m ri 

- ,/ 2 E E CT^ m (h) / WK- 
2ttV1-cos 2 p £=2imn J c(r;P ) 

Applying 



y, n (r')dr' = P,(r)F/(r) 
(Freeden et al., 1998, p. 64) with r = cosp finally yields 

(W/)(h,r;p) = E^(cosp)Ecr^ m (h)^(r), (30) 

£=2i m,n 

(cf. Bunge, 1969, p. 45; Bunge, 1982, p. 74). 

In case of an orientation probability function /(arccos(Sc(gQg))) 
which is radially symmtric with respect to a given qo, Eq. (|3T)|) simplifies 
further to 

(W/)(h,r;p) = E^(cosp)ECT^ m (h)^(r) 

= V^7^C^(cosp)P/(cos77) (31) 

as the Radon transform becomes a function of 77 = arccos(gohgo ■ r ) = 
d(q , C) only, where C denotes the great circle representing all rotations 
mapping h on r. It should be noted that p and 77 commute, i.e. 



(Wf)(h,r;p)\ qohq *. r=cosr] = (Wf)(h,r;r))\ qohg * 



r=cos p 



Integrating the radially symmetric spherical Radon transform over a 
small circle with angle p with respect to r with (/oh^o ' r = cosr 7 is 
equal to integrating the spherical Radon transform over a small circle 
with angle 77 with respect to r with qohq^ ■ r = cos p. In terms of the 
generalized spherical Radon transform we find 
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Proposition 7. For a radially symmetric function f 

(TZ^ P ^ /)(Ch, r )|g hg5-r=cos?j = (7?.^ f) (Ch,r) \q hq%T=cosp 

Sampling (W/)(h, r;p) for discrete values of h, r G S 2 and p G 
(0, 7r/2) gives rise to data Wi = w(hi, p^) and a system of linear equa- 
tions 

Wi = Y, Pi(cos Pi ) CrYr(^i)Y e n (ri), i = 1, . . . , I 

1=21 m,n 

which may allow a solution for the C^-coefficients with even £. 

4.3 Convolution kernels and radial basis functions 

Assume that the orientation probability density function shall be mod- 
eled by the superposition of "components" or kernels Ki(q), i.e. 

f(q) = ^ Kjfa ff) = E X J K fa> «)■ 
j j 

Usually in practice, we choose radially symmetric kernels which give rise 
to 

/(?) = ^A iJ fsT(arccos(Sc(g*g))). 

3 

Then the corresponding generalized spherical Radon transform Wf 
is the superposition of the correspondingly generalized spherically Radon 
transformed kernels kj = WKj 

(W/)(h,r;p) = ^A i (WA i )(h ) r;p) 

i 

= ^A^Qi^p), 

which may be fitted to the experimental data 

Wi ~ ^ A* fc(hj, r 4 ; p<) 

i 

in some sense, e.g. in the sense of a Hilbert-Sobolev norm as developed 
in (Boogaart et al., 2005). Then, an approximate orientation probability 
density function explaining the data is given by 



/(g) = EA*K(arccos(Sc(g;g))). 
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5 Examples 



In the following we provide some formulae for the Abel-Poisson (in prob- 
ability: Cauchy) and the de la Vallee Poussin kernel, their Radon trans- 
form and their twofold Radon transform. 

In texture analysis the Abel-Poisson kernel is referred to as Lorentz 
standard function (Matthies et al., 1987, p. 98; Matthies et al., 1990, 
p. 477), and the formulae were actually taken from there. Obviously, 
they were initially not related to reproducing kernels and their twofold 
spherical Radon transform. The de la Vallee Poussin kernel has been 
introduced into texture analysis because of its harmonic series expansion 
is finite (Schaeben, 1997; 1999). 

For the symmetrical kernel K(pi,p 2 ) = K(plp 2 ) defined on § 3 x S 3 , 
the variable uj = 2 arccos(Sc(p^ 2 )) denotes the angle of the rotation of 
p\p 2 ] for the Radon transformed kernel TZK&^p, h, r) defined on S 3 x § 2 x 
S 2 , the variable rj denotes the angle Z(php*, r). The two variables rji and 
r}2 of the twofold Radon transformed kernel 1Z1ZK(h.i, ri, h 2 , r 2 ) defined 
on (8 2 x § 2 ) x (§ 2 x § 2 ) correspond to the angles Zhih 2 and Zrir 2 . 
The Gegenbauer respectively the Legendre coefficients of the kernels are 
denoted ag. 

In the table we have used the following notations for special functions: 
B Beta function, 2 Fx hypergeometric function, and T Gamma function,. 
For the Abel-Poisson kernel we have 



K 
TZK 



a e = (2£+1)k z 

1-K 2 1-K 2 



2 L(l - 2/tcos(o;/2) + k 2 )2 (1 + 2kcos(cj/2) + k 2 )2 
1-k 4 



1 -2k 2 COS77 + K 4 ) 3 / 2 



2 1-K 2 2D . 

1Z1ZK = E - — - 

tt(C-D)VcTD y C + D> 

where C = 1 — 2k cos 771 cos r] 2 + k 2 and D = 2k sin 7]\ sin ?? 2 
Analogously, for the de la Vallee Poussin kernel 

at = (2B(|«+i))- 1 [ft(*)-5' m («)] 
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where S £ (k) = ££ =0 (-l) fe Ql) B ( k + h K + i ~ k + D> and further 

B(3/2,« + l/2) 1 7 ; 
ftA = (1 + k) cos(t?/2) 2k 



7T cos r/i cos 772 T(| + k) 



where A = 1 + cos(?7i + 772) and 5=1 + cos(?7i — 772). 



6 Conclusions 



The essential role of the probability density of the angle distribution 
for the inverse spherical "totally geodesic" Radon transform has been 
clarified by purely geometric arguments. It is identified with the gen- 
eralized spherical Radon transform which in turn is identified with the 
"spherically translated" spherical Radon transform. Of particular in- 
terest is that the twofold spherical Radon transform of a symmetrical 
kernel function is again its corresponding angle probability density func- 
tion. Practical methods of inversion in terms of harmonics or radially 
symmetric basis functions are sketched. Thus, our contribution is also a 
tribute to the late Hans- Joachim Bunge (1929 - 2004), who introduced 
the angle distribution into "quantitative texure analysis". The problem 
whether the inversion of the generalized spherical Radon transform is 
better conditioned than the inversion of the spherical radon transform 
is postponed to a future contribution as it requires a detailed analysis 
of the experiment to collect integral radiation intensity data. 
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